Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN MfeK style file vl.4) 



The correlation function of radio sources 



A. J. Loan, 1 J. V. Wall 2 and O. Lahav, 1 

1 Institute of Astronomy, Madingley Road, Cambridge, CBS OH A 

2 Royal Greenwich Observatory, Madingley Road, Cambridge, CB3 0EZ 



Os 
Os 

o 

Q 

OS 



> 
o 
os 

(N 
OS 

a: 

6 

c3 



Received ; accepted , 1996 



ABSTRACT 

We investigate the large-scale clustering of radio sources in the Green Bank and Parkes- 
MIT-NRAO 4.85 GHz surveys by measuring the angular two-point correlation function 
w(9). Excluding contaminated areas, the two surveys together cover 70 per cent of the 
whole sky. We find both surveys to be reasonably complete above 50mJy. On the 
basis of previous studies, the radio sources are galaxies and radio-loud quasars lying 
at redshifts up to z ~ 4, with a median redshift z ~ 1. This provides the opportunity to 
probe large-scale structures in a volume far larger than that within the reach of present 
optical and infrared surveys. We detect a clustering signal w(8) w 0.01 for 9 — 1°. By 
assuming an evolving power-law spatial correlation function in comoving coordinates 
£(r c , z) — (r c /ro)~ 7 (1 + z) J ~( 3+e \ where 7 = 1.8, and the redshift distribution N(z) 
of the radio galaxies, we constrain the rg-e parameter space. For 'stable clustering' 
(e = 0), we find the correlation length tq m 18ft. _1 Mpc, larger than the value for 
nearby normal galaxies and comparable to the cluster-cluster correlation length. 

Key words: galaxies: clustering - radio galaxies - large-scale structure 



1 INTRODUCTION 



Surveys of galaxies selected by their optical or infrared emis- 
sion reveal the rich structure of the universe. This clustering 
structure can be quantified by using statistical techniques, 
such as the two-point correlation function in two (e.g. APM: 
Maddox et al. 1990) or three (e.g. IRAS: Fisher et al. 1993) 
dimensions. However, these surveys only probe the local uni- 
verse: for example, the median redshift of the APM survey 
is z ~ 0.1. In contrast, radio galaxies and radio-loud quasars 
can be detected over significant cosmological distances, up 
to redshifts z ~ 4. However, such bright sources of radio 
emission represent a small fraction of all galaxies, sampling 
the large volume within which they can be observed much 
more sparsely than normal galaxies. 

The prevailing consensus has long held that the distri- 
bution of radio galaxies is almost isotropic, with any struc- 
tures barely detectable. Webster (1976) studied the dis- 
tribution of over 7500 sources from 3 separate radio sur- 
veys at 148 MHz, 408 MHz and 1.4 GHz, concluding that the 
data were consistent with a uniform random distribution to 
within a few per cent. Shaver & Pierre (1989) demonstrated 
that the distribution of nearby (2 < 0.02) radio galaxies was 
flattened towards the Supergalactic Plane. More recently, 
Peacock & Nicholson (1991) analyzed a redshift survey of 
300 nearby (z < 0.1) radio galaxies, measuring a power- 
law spatial correlation function with a correlation length 
ro = 11.0±1.2 ft -1 Mpc, much larger than that for optically- 



selected galaxies (ro ~ 5/i _1 Mpc - Peebles 1980). Benn 
& Wall (1995) showed how the apparent isotropy of radio 
sources may be used to estimate the largest scales of cellu- 
lar structure; further background and references are given 
therein. 

The clustering of radio sources in the 4.85 GHz Green 
Bank survey of the northern hemisphere (87GB: Gregory 
& Condon 1991) has been studied recently, making use of 
the angular two-point correlation function w(8) (Kooiman, 
Burns & Klypin 1995; Sicotte 1995). Both of these studies 
find that there is clearly detectable clustering signal in the 
87GB survey. A preliminary correlation function analysis 
of the first ~15 per cent of the FIRST survey also finds 
significant clustering of radio sources (Cress et al. 1996). 

This paper presents a new study of the large-scale struc- 
tures traced by the radio sources in the 87GB survey. We also 
present a correlation function analysis of the Parkes-MIT- 
NRAO survey (PMN: Griffith & Wright 1993). The PMN 
survey used much the same method as the 87GB survey to 
produce a complementary 4.85 GHz survey covering most 
of the southern sky. With appropriate selection of angular 
coverage and flux-density limit, the two surveys prove to 
be quite well matched, covering most of the celestial sphere. 
Nevertheless, we analyze the surveys separately to avoid any 
problems arising from small systematic differences in source 
detection procedure and flux-density limit. 

In Section ^ we examine the 87GB and PMN surveys to 
demonstrate their suitability for correlation function analy- 
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Figure 1. The 87GB (northern hemisphere) and PMN (southern hemisphere) source catalogues in Equatorial coordinates, Aitoff 
projection. The Galactic Plane (|6[ < 10°) and areas of poor coverage or confusion (see Table [l]) have been excised. The varying source 
density reflects the declination dependence of the flux limits of 87GB and PMN, quantified in Figure 0. 



sis. In Section |3j we present the results of our analysis, while 
Section ^ outlines the formalism that we use to estimate the 
parameters of the correlation function in three dimensions. 
In Section ^ we discuss the implications of our results and 
compare them to those of other studies. 



2 THE RADIO CATALOGUES 

2.1 The Green Bank survey (87GB) 

The 87GB source catalogue was compiled by Gregory & 
Condon (1991) from observations taken in 1987 October 
by Condon, Broderick & Seielstad (1989) using a seven- 
beam receiver attached to the (now defunct) NRAO 300 foot 
(91 m) radio telescope in Green Bank, West Virginia, USA. 
The telescope was tracked back and forth along the local 
meridian at ±10° min -1 , the diurnal rotation of the earth 
eventually providing a Nyquist spacing of a half beam- width 
between adjacent observations. The data were reduced auto- 
matically using a 5-<r source-detection criterion to produce 
a catalogue of 54,579 sources covering the declination band 
(0° < 8 < +75°) for all right ascensions, except for four 
small confused or noisy areas. The 87GB survey has a total 
sky coverage of ~ 6.0 steradians. 



2.2 The Parkes-MIT-NRAO survey (PMN) 

The PMN source catalogue was compiled from observations 
taken in 1990 June and 1990 November with the ANRAO 
64 m radio telescope in Parkes, NSW, Australia, covering 
the declination range —87.5° < S < +10.0° for all right 
ascensions. These observations followed the same observing 



Table 1. Areas from the 87GB and PMN surveys included in 
and excluded from our analysis. 

Included: 
87GB: 
PMN: 



Excluded: 
Galactic Plane: 
Solar interference: 
Extended sources: 



75° > 8 > 1° 
0° > 8 > -28° 
-38° > 8 > -85° 

|6| < 10° 

6° > 8 > 1° and \a - 205° < 20° 
8 < -60° and \a - 80° | < 10° 
8 > -10° and \a - 205° | < 30° 
-30° > 8 > -50° and \a - 200° < 10° 



procedure using the same seven-beam receiver, and were re- 
duced using an algorithm similar to that used with 87GB. 
However, a slightly different source-detection criterion was 
adopted: 90 per cent reliability - equivalent to ~4.4-a - 
rather than a 5-er. 

The source catalogue was split into 4 zones according 
to the elevation of the telescope (see Table 1 of Wright et al. 
1994): Southern, Tropical, Equatorial and Zenith (respec- 
tively: Wright et al. 1994; Griffith et al. 1994; Griffith et 
al. 1995; Wright et al. 1996). However, because the Parkes 
telescope tracks very slowly in the 10° wide region near the 
local zenith, the data from the Zenith zone were observed 
using a different strategy to the other three zones and is of 
a much lower quality. For these reasons, the Zenith zone is 
excluded from our analysis, leaving the other three zones 
which together cover an total area of ~ 6.4 steradians. 
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Figure 2. Histogram of the surface density of radio sources in 
the 87GB (8 > 1°) and PMN (S < 0°) source catalogues as a 
function of declination, for a range of flux-density limits (35, 50, 
70 and 100 mjy). 



2.3 Comparing the 87GB and PMN surveys 

Before the 87GB and PMN source catalogues can be used for 
correlation function analysis, we need to show that they are 
well matched, with no large gradients caused by instrumen- 
tal or other systematic effects. We note that real large-scale 
structures, such as the Supergalactic Plane, will introduce 
inhomogeneities in these catalogues on large angular scales. 
We will address this issue elsewhere (Baleisis et al. 1996) but 
here assume that real large-scale structures appear on small 
angular scales and use whole-sky uniformity as our criterion 
for completeness. 

The 87GB and PMN surveys overlap in the declination 
band 0° < 6 < 10° (» 0.59 steradians). Griffith et al. (1995) 
examined this overlap region, finding agreement within the 
quoted errors between both the published flux densities and 
the positions of the sources. However, both Griffith et al. 
(1995) and Baleisis (1996) report a worrying 2 per cent dif- 
ference between the flux calibration of the 87GB and PMN 
catalogues, in the sense that sources in the PMN survey 
are assigned a systematically lower flux density than those 
in the 87GB survey. For this reason, and to highlight any 
discrepant results, we treat the surveys separately. We note 
that such a small calibration offset will not significantly af- 
fect the results of our analysis on small scales, as we discuss 
below. 

The flux-density limit of the 87GB survey varies with 
declination, increasing from ~25 mjy at high declinations 
(6 > 60°) to ~40mJy at the equator (Gregory & Condon 
1991), for two main reasons. First, there is an elevation de- 
pendence of the sensitivity of the 300 foot Green Bank tele- 
scope, and secondly, adjacent observing tracks become in- 
creasingly further apart at lower declinations. The Parkes 



telescope does not suffer from the problem of decreased sen- 
sitivity at low elevations. However, the increased coverage 
due to adjacent scans lying closer together at far-southern 
declinations (S < —70°) causes the flux-density limit to de- 
crease from ~50mJy near the local zenith (5 ~ —33°) to 
~40 mjy near the southern equatorial pole. 

We use 87GB in the northern hemisphere (5 > 1°) 
because 87GB has a flux limit which is lower than that 
for PMN over the region of overlap. We discard the region 
within 10° of the Galactic Plane to avoid contamination by 
Galactic radio sources, and we also discard the area where 
there is solar interference in the 87GB survey, together with 
three areas with confusion from extended sources in the 
PMN catalogue. The discarded areas are shown in Table ^. 

The remaining sources, shown in Figure [j], form the cat- 
alogue which we use in our correlation analysis. Figure |l| still 
shows a strong variation in source density with declination 
due to changes in the local flux-density limit. To counter- 
act this effect, we must impose a higher flux limit over both 
surveys. 

We plot the surface density of radio sources in the 87GB 
and PMN catalogues as a function of declination in Figure ^ 
for flux-density limits of 35, 50, 70 and 100 mjy in order 
to find a safe lower flux-density limit, which will provide a 
source catalogue with uniform sensitivity across the whole 
area. Both surveys are largely incomplete at the 35 mjy flux 
limits, with large density gradients caused by the sensitivity 
problems discussed above. The surface density of sources 
becomes increasingly more uniform as the limit is raised 
from 50 mjy, through 70 mjy to 100 mjy. Above 50 mjy, the 
surveys appear to be sufficiently evenly matched over the 
whole sky to allow structure on small angular scales to be 
studied using correlation techniques. We will also calculate 
correlation functions at a limit below 50 mjy to demonstrate 
the effect of variable sensitivity on the measured clustering. 

There remains some suggestion of larger-scale structure 
present in the combined 87GB/PMN source catalogue which 
we examine elsewhere (Baleisis et al. 1996). 



3 THE ANGULAR CORRELATION 
FUNCTION 

The two-point correlation function is a measure of how clus- 
tered the sources are compared to a random Poisson dis- 
tribution (e.g. Peebles 1980). The angular two-point corre- 
lation function w(9) gives the excess probability of finding 
two sources in the solid angles <5fii and 6Q2 separated by 
an angle 6. One way to estimate this function (Hamilton 
1993) is to compare the distribution of real galaxies to the 
distribution of points in a random catalogue with the same 
boundaries, that is: 



w(9) 



DD.RR 
{DRY 



(1) 



where DD, RR and DR are the numbers of data-data, 
random-random and data-random pairs separated by the 
distance 9 + 89. The estimation of RR and DR requires a 
catalogue of 'objects' scattered uniformly over an area with 
the same selection properties as the data catalogue (i.e. with 
the same angular boundaries). The number of 'objects' in 
the random catalogue may be very much greater than that 
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Figure 3. Angular correlation function w{&) for the 87GB survey, with flux-density limits of 100 mJy, 70mjy, 50mjy and 35mjy. In 
each case the dashed line shows the data in 0.2° bins, and the solid line in 1° bins. Poisson error bars are shown on the line for 0.2° bins. 
The excess positive signal at small angles (9 < 2°) increases in amplitude at higher flux-density limits. 



in the real catalogues to reduce the errors in the estimates 
of DR and RR. 

We have calculated the angular two-point correlation 
functions for the 87GB and PMN surveys separately, to 
highlight any discrepant results and also to avoid problems 
if the catalogues prove to be mismatches. The calculated 
correlation functions for the 87GB and PMN catalogues are 
shown in Figure |^ and Figure ^, respectively. These figures 
show an excess positive signal (i.e. w(6) > 0) at small an- 
gles (6 < 2°) with increasing amplitude at higher flux lim- 
its. This increased amplitude arises from the 'washing out' 



of structures by foreground or background sources in sur- 
veys with a lower flux-density limit. A shallower survey is 
expected to show a higher amplitude of the correlation func- 
tion since this 'washing out' is less severe. However, a sur- 
vey with a higher flux-density limits contains fewer sources, 
leading to larger errors in the estimate of w(0) which may 
conceal the greater expected amplitude. 

The correlation function for the 35 mjy-limit PMN cat- 
alogue, fourth panel of Figure ^, shows the effect of large- 
scale density gradients in the survey noted from Figure g. 
These gradients are caused by choosing sources at a flux- 
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Figure 4. Angular correlation function w{&) for the PMN survey, for flux-density limits of 100 mJy, 70 mJy, 50 mJy and 35 mJy. In each 
case the dashed line represents 0.2° bins and the solid line 1° bins. Poisson error bars are shown for the 0.2° binning. There is excess 
positive signal at small angles (9 < 2°); the 35mjy plot shows the effect of the large density gradients below a flux-density completeness 
limit. 



density limit below that of the catalogue completeness, pro- 
ducing a large error in the measured correlation function. 

We assume the usual power-law form for the angular 
correlation function, such that 

w{6) = AO 1 -"'. (2) 

If we assume 7 = 1.8, as found for local galaxies (Peebles 
1980), then only the amplitude A of w(0) remains to be 
estimated from our measurements. We have estimated this 



amplitude by a simple least-squares procedure, by summing 
over the i bins of our calculated correlation function: 

A t _ E^gr , 3) 

where Wi is the observed value of w(0) in the bin at angular 
separation 8i. Note that different bins of the observed w{9) 
are correlated. Figure ^| shows our estimates for the ampli- 
tude of the correlation function for various binning schemes 
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Figure 5. Estimated correlation function amplitude A in the ex- 
pression w(8) = AO 1 ^ 7 for a variety of different binnings, for each 
chosen flux-density limit. 87GB is shown with a dashed line, PMN 
with a solid line. The lowest flux-density limit, 35 mjy, must be 
regarded as suspect due to incompleteness. The remaining suggest 
a systematic trend with higher amplitude at higher flux limit. 



and flux limits. Our estimates indicate that the amplitude 
A has a value between 0.005 and 0.015. 



4 THE SPATIAL CORRELATION FUNCTION 

4.1 The cosmological Limber's equation 

Limber's equation (Limber 1953) relates the spatial corre- 
lation function £(r, z) (in 3 dimensions) to the angular cor- 
relation function w{6) (in 2 dimensions), given the selection 
function of the sample and a cosmological model. The cos- 
mological Limber's equation can be written as (e.g. Section 
56 in Peebles 1980; Baugh & Efstathiou 1993) : 



u,(0) 



2 J °° J °° x A F- 2 4>\x)ti{r, z)dxdu 
[J™ x 2 F-icj>(x)dx\ 2 



(4) 



where a; is a comoving coordinate and F takes into account 
the different possible world geometries (e.g. for Qq = 1, Ao = 
0, then F = 1). The selection function 4>{x) satisfies the 
following relations for the mean surface density in a survey 
of solid angle f2 3 : 

Af= x 2 (f>(x)dx = — / N(z)dz, (5) 
Jo fis Jo 

where N(z)dz is the number of galaxies in the survey in the 
redshift shell (z, z + dz). 

Hereafter we assume a flat universe with Qo = 1, F = 1, 
for which the comoving distance is 



x=%-[l-(l+z) 
Ho 



-1/2, 



(6) 



where Ho = 100 h~ J kmsec - Mpc -1 is the Hubble ex- 
pansion parameter. We also assume a power-law redshift- 
dependent 3-dimensional correlation function: 



£(r,«) = 



(0 



(1 + *) 



(7) 



where r is a proper coordinate, 7~o is the correlation scale 
length at redshift z — 0, and e describes the redshift evolu- 
tion of the spatial correlation function. In comoving coordi- 
nates, this can be written as: 



l + z 



,7-(3+e) 



(8) 



where r c = r(l + z) is a comoving coordinate. Specific val- 
ues of e can be interpreted as follows: e — corresponds to 
constant clustering in proper coordinates, i.e. 'stable clus- 
tering'; e = 7 — 3 = —1.2 implies constant clustering in 
comoving coordinates; e = 7— 1 = 0.8 implies growth of 
clustering under linear theory (e.g. Peebles 1980; Treyer & 
Lahav 1996). 

For our purpose it is more convenient to express Lim- 
ber's equation in term of N(z). Using Equation ^ and Equa- 
tion ^| we can write: 



! (l + z 



,+3/2 



(9) 



Substituting this into Equation ^| and assuming the small 



angle approximation r 
that: 

(as ( r o H o 



(u 2 /F 2 + x 2 9 2 ) 1/2 /(l + z) we find 



(10) 



[l + z) 



1-7 



n 2 {z) + 



[J?dzN(z)Y 



where H~, 



• - ^ {\) T [^T-) /r(^), with F the Gamma func- 
tion. Assuming 7 = 1.8, Hi.g = 3.68. 

We emphasize that these equations are only valid for 
a power-law correlation function of the form of Equation ^ 
and fio = 1 universe. A more general expression is given in 
terms of a double integral and Fourier transforms by Baugh 
& Efstathiou (1993) and Treyer & Lahav (1996). 

Since 7 is assumed a priori to be 1.8, the same as that 
for normal galaxies, and N(z) is determined from studies 
of the radio-source population, the only free parameters to 
be deduced by comparison with the observed w(ff) are the 
correlation scale length ro and the clustering evolution index 
e. Note that the amplitude of N(z) does not affect Limber's 
equation; only the functional form of N(z) is important. 



4.2 The redshift distribution N(z) 

The prediction of the angular two-point correlation func- 
tion w(8) from the spatial correlation function £(r) requires 
knowledge of the redshift distribution of the sources N(z) 
complete to the designated flux-density limit. No complete 
sample of radio sources has been fully identified at the flux- 
density level of 50mJy. Despite this, Dunlop & Peacock 
(1990) have shown how a synthesis of the statistics of ra- 
dio sources (including source-counts together with incom- 
plete identifications and measured redshifts) may be used 
to determine the distribution of radio sources with redshift, 



© 0000 RAS, MNRAS 000, 000-000 



Correlation function of radio sources 7 




40 



Redshift 

Figure 6. Redshift distributions N(z) of the radio-source pop- 
ulation selected at 4.85 GHz, shown for each flux-density limit at 
which w(9) has been calculated. The solid curves represent the 
average of six models taken from Dunlop & Peacock (1990); see 
text. The dotted curves for each limit represent the 6 models 
from which the average was formed, indicating the uncertainty at 
each flux-density level due to incomplete or statistically-limited 
redshift data. 



from which N(z) may be derived at any flux density. In their 
analysis it was assumed that there are two independent pop- 
ulations of radio sources, namely the flat-(radio)spectrum 
and steep-spectrum objects. The former are predominantly 
identified with QSOs, the latter with radio galaxies; a more 
detailed overview of the radio-source populations is given by 
Wall (1994). Treating the data for flat- and steep-spectrum 
sources separately, Dunlop & Peacock fit polynomial ex- 
pansions to describe the luminosity function and its epoch- 
dependence for each population. They did so from seven dif- 
ferent starting points, the scatter amongst the seven models 
providing indication of the statistical definition of the space 
density at any point over the P — z plane. 

From these models we have calculated the N(z) at 35, 
50, 70 and 100 mjy, summing the space densities for the 
flat-and steep-spectrum populations. At each flux-density 
limit, we then take the mean from the various models as the 
the canonical value for N(z). (In this procedure we omitted 
model 5, which has a problem at the lower flux densities 
in that a dominating spike appears in the space densities 
and N(z) close to the redshift cutoff which is specified as a 
starting condition for this particular formulation.) The N(z) 
produced by this procedure is shown in Figure |^. This figure 
indicates that the median redshift for the radio sources in 
the 87GB/PMN survey is at z ~ 1. Figure ^| also shows 
that the average redshift decreases as the lower flux limit is 
increased. 

We note that a 2 per cent flux mismatch between the 
catalogues has been reported (Griffith et al. 1995; Baleisis 
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Figure 7. Predicted loci of constant A in the correlation length 
ro and clustering evolution index e plane, for S4.85 > 100 mjy 
(solid lines) and S4.85 > 50 mjy (dotted lines). The values of ro 
and e which may be deduced from a single measurement of A are 
strongly correlated. 



1996). This mismatch would mean that the flux limits given 
for the correlation functions estimated above are slightly in- 
correct, and that we should use a different N(z) correspond- 
ing to the true flux in our Limber's inversion. However, as 
we can see from Figure ^| the shape of N(z) varies very lit- 
tle from one flux limit to another, and so there will be little 
effect on the small angular scales concerned. 



4.3 Estimating the spatial correlation function 

t(r,z) 

Armed with our estimate for the amplitude of the angular 
correlation function from Section H (A = 0.010 ±0.005) and 
our estimate of N(z) from Section }4.2| , we can now compare 
this observed estimate to the predictions from our model of 
the radio galaxy population for a range of parameters, via 
equations 2 and 10. Figure |^ shows lines of constant w(6) 
amplitude as a function of the spatial correlation length ro 
and the evolution index e at a flux-density limits of 50 mjy 
and 100 mjy. There is a strong correlation between the val- 
ues of ro and e which can be inferred from a single mea- 
surement of A. Assuming a value for one parameter then 
provides an estimate for the other parameter. For example, 
consider w(8 = 1°) = A = 0.01 for a flux limit of 50mJy. If 
we also assume that clustering is constant in comoving coor- 
dinates, i.e. e = 7 — 3 = —1-2, then ro = 13 h~ x Mpc, while 
if clustering is 'stable', i.e. e = 0, then ro = 18 h~ Mpc. The 
range of acceptable values of A (0.005 to 0.015) translates 
into a uncertainty of a few Mpc in ro. 



© 0000 RAS, MNRAS 000, 000-000 



8 A.J. Loan, J. V. Wall and O. Lahav 



5 DISCUSSION 

Using two large-area 4.85 GHz radio surveys with median 
redshift z ~ 1, we have detected a correlation signal w(8) « 
0.01 at angles 6 ~ 1° for radio sources brighter than 50mJy. 
Adopting the best estimate of redshift distribution for the 
sources, we find the present-day correlation length to be in 
the range 13 < ro < 18 h~ Mpc for evolution index in the 
range —1.2 < e < 0, with uncertainty in the measurements 
permitting errors as large as ~ 5/i _1 Mpc. The values of 
ro and e inferred from a measurement of w(9) are strongly 
correlated. 

Previous studies by Kooiman et al. (1995) and Sicotte 
(1995) measured angular correlation functions which were 
largely in agreement with each other and with our results. 
Kooiman et al. and Sicotte both attempted to use their re- 
sults to constrain the spatial correlation function of radio 
galaxies, but reached contradictory results. Sicotte found 
good agreement between his model and his measured cor- 
relation function, with a best fit value of e = —1 and 
ro = 11 h^ 1 Mpc. This result is in good agreement our re- 
sults, even though Sicotte used an earlier (Condon 1984) and 
significantly different estimate of N(z). (We repeated our 
analysis using the Condon N(z) and obtained remarkably 
similar results; there is little sensitivity to the form of N(z).) 
This is also in accord with the value of ro = 11 h^ 1 Mpc 
measured by Peacock & Nicholson (1991) from low redshift 
(z < 0.1) radio galaxies. 

On the other hand, Kooiman et al. found that the am- 
plitude of their model was 70 times smaller than the mea- 
sured correlation function, assuming 7 = 1.8, e = and 
ro = ll/i _1 Mpc. This is clearly inconsistent with our re- 
sults, and also with those of Sicotte. 

Taken at face value, for a realistic range of the evolution 
index e, we find a present-day correlation length ro which 
is far larger than that for normal optical galaxies (ro ~ 
5 ft" 1 Mpc, Peebles 1980), IRAS galaxies (r ~ 4/i" 1 Mpc, 
Saunders, Rowan- Robinson & Lawrence 1992), and QSOs 
(r ~ 6 h- 1 Mpc for e = -1.2, Shanks & Boyle 1994). How- 
ever, our values of the correlation length for radio galaxies is 
comparable to that for clusters of galaxies (ro ~ 20 h^ 1 Mpc, 
Bahcall 1988). This may reflect the preference of radio galax- 
ies at z > 0.5 to reside in high-density environments (Bah- 
call & Chokshi 1992; Yee & Ellingson 1993 and references 
therein). 

The uncertainty in our measurement of the correlation 
function could be reduced by measuring the redshifts for at 
least a random subset of the catalogues and cross-correlating 
this subset with the rest of the sample (e.g. Saunders et al. 
1992). These redshifts would also improve our knowledge 
of the redshift distribution of the radio galaxy population 
N(z), and indeed of the nature of the population (QSO - 
radio-galaxy - BL Lac object) at this flux-density level. 

Our analysis demonstrates that large-scale structure 
studies based on new radio surveys nearing completion will 
provide decisive estimates of distant structure and its evolu- 
tion; and we have provided the appropriate formalism. On- 
going surveys with the VLA - FIRST (Becker, White & 
Helfand 1995) and NVSS (Condon et al. 1996) - will yield 
of order 10 sources over the sky, thus overcoming the lim- 
itations imposed by counting errors in the estimated pair- 
counts from the current 4.85 GHz surveys. The first corre- 



lation analysis of a subset of the FIRST survey (Cress et 
al. 1996) indicates that the full survey will provide strong 
constraints on our models of the radio source population. 
Furthermore, positional accuracy of these new surveys will 
be such that direct spectroscopy at the catalogue radio po- 
sitions can be used to obtain redshifts. Knowledge of the 
redshifts of even a random subset of these catalogues will 
allow the use of cross-correlation techniques to improve the 
estimate of the correlation function in three-dimensions (e.g. 
Saunders et al. 1992). 
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